Model for vacancy-induced dsup Oferromagnetism in oxide compounds. 
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We propose a model with few parameters for vacancy-induced ferromagnetism based on a corre- 
lated model for oxygen orbitals with random potentials representing cation vacancies. For certain 
potentials, moments appear on oxygen sites near defects. Treating the randomness exactly, we 
calculate the magnetic couplings between moments, the Curie temperature and spin and charge 
densities as a function of the potential, the density of vacancies, and correlation strength. For 
physically reasonable parameters this predicts Curie temperatures well above room temperature for 
small concentrations of vacancies. We discuss our results in relation to questions of stability and 
reproducibility raised in experiments. To circumvent the difficulties of controlling intrinsic defects, 
we propose specific non-magnetic host doping that could be, for example, substituted for cations in 
Hf0 2 or Zr0 2 . 



present letter is to provide a model for vacancy induced 
ferromagnetism explaining Curie temperatures at room 
temperatures. This can also help understand instabili- 
ties, and propose new materials where controlled substi- 
tutions could replace the vacancies. Point defects have 
been suggested several times as a mechanism for local- 
ized moments 0, 0, Hij l and ferromagnetism in oxides. 
While oxygen vacancies were originally suggested as the 
source of magnetism in HfC>2 @, density functional cal- 
culations including various defects found no evidence 
for moments around such defects. A large magnetic mo- 
ment was, however, found in partially depleted oxygen 
orbitals around Hf vacancies. Couplings were estimated 
between two defects confined within a supercell. This 
is not sufficient to estimate properly the Curie tempera- 
ture as a function of vacancy concentration which require 
couplings at all distances. We propose a theory where de- 
pendence on concentration and also the effective doping 
of different defects can be studied. We start from the idea 
of vacancy- induced moments and construct a model, with 
only few free parameters, in which we can estimate not 
only the temperature dependence of local moments, but 
also the couplings required to maintain long-range ferro- 
magnetic order at the low concentrations relevant to the 
experiments. We can also calculate the thermodynam- 
ics and dynamical spin correlations for future systematic 
experimental investigations. We formulate the problem 
for the cases of cation vacancies in oxides such as Hf02, 
ZrC>2 and CaO, but the aim is to explore what may be 
a more geneal phenomenon, in which a vacancy or sub- 
stitutional defect creates an extended magnetic moment 
on neighboring atomsWhile there are similarities to the 
theory of magnetism in diluted magnetic semiconductors 
but there the substitutional impurities immediately pro- 
vide localized magnetic moments which interact via the 
itinerant carriers of the (doped) host. Here || the mo- 
ment is induced on several non-magnetic oxygen atoms of 



A number of exciting, but puzzling, effects when ferro- 
magnetism is seen in unexpected places, are apparently 
unrelated to traditional transition metal magnetism: in 
oxides such as HfC-2, ZrC-2 lj, CaO, ZnO, and related ma- 
terials such as hexaborides[2 (CaBg) and even irradiated 
graphite HQ. Thin films of materials which in bulk have 
neither magnetic moments nor magnetic order, may be 
ferromagnetic well above room temperature^, [g. Hf02 
or Zr02, are both wide-band insulating material with 
high dielectric constant; the possibility of making them 
ferromagnetic could widen their possible application in 
the field of spintronics. There are analogies to the dis- 
covery that the doping by a small amount of mag- 
netic impurities could lead to relatively high Curie tem- 
perature in III-V semiconductors. There are, apparently, 
however too few magnetic impurities in the Zr02 films, 
for example, to explain the magnetism and the cation 
(nominally Zr 4+ ) should be in a non-magnetic d° con- 
figuration. The term " d°" ferromagnetism was coined to 
describe this general phenomenon. There are suggestions 
that the ferromagnetism is related to intrinsic cation or 
anion vacancies, and that this may be a path to new 
ferromagnets|], Ej-Band structure calculations [lfl lrH | 
supported the idea that vacancies can induce local mo- 
ments on neighboring atoms. What is lacking, espe- 
cially as the experimental results are still unstable, is 
a quantitative theory of long-range ferromagnetic order 
of any such moments. This can define the important 
parameters, for example the nature and concentrations 
of vacancies and their effective doping. The aim of the 
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the original host and their is no clear separation between 
the holes forming the moment and those mediating the 
interactions. 

We describe this via a single correlated band of oxygen 
orbitals, with random potentials representing the influ- 
ence of the randomly introduced vacancies on neighbor- 
ing cation orbitals. The oxygen atoms are on a regu- 
lar lattice, taken to be cubic for simplicity. In the pure 
host, the parameters can be those calculated for the filled 
oxygen-dominated bands. The cations in the pure mate- 
rial are assumed to be at the center of each elementary 
cube of the lattice. A concentration x of cation vacan- 
cies is modeled by choosing randomly the positions of 
cube centers. Thus a "defect" does not alter the regular 
geometry of the oxygen lattice but introduces attractive 
potentials on the 8 neighboring oxygen orbitals of the 
cube center chosen. The usual Hubbard Hamiltonian on 
a regular lattice has an extra term which is a random, 
but locally correlated, potential: 

wo = -5>.*' (4 a <' + < -4° 

+ ^t/nJnJ--K(nJ+nj-) 

i 

c^ (c) operators create (destroy) electrons in the oxygen 
orbitals. Vi are attractive on sites adjacent to the va- 
cancies, corresponding to attraction of the holes. Near- 
est neighbor hopping tij = t everywhere except between 
pairs of sites neighboring the same vacancy for which 
tj j = t' . Defects in Hf/ZrC>2 have two inequivalent sets 
of neighboring oxygens (01 and 02) [13. To ease compar- 
ison with Local Density Approximation density of states 
we took Vi = V on four neighbors and Vi = y on the 
other four. V is a free parameter. The hole density is 
taken to be a variable nt ■ In the simplest picture (only 
cation vacancies) rih is Ax for HTO2 and Zr02 (cation 
4+) and 2x for CaO or ZnO (cation 2+). In numerics we 
quote U and V in terms of the host bandwidth W = 12t. 
For Zr02 we take W — 7 eV from the calculated band- 
width of the p band just below the Fermi level 01 an d 
U ~ 3 eV as estimated for the Hubba rd p arameter for 
oxygen p orbitals U pp in copper oxides [14|. 

We treat the correlated disorder exactly as the unre- 
stricted Hartree-Fock (UHF) approximation is applied in 
real space for each configuration of disorder, sampled over 
many configurations. Comparison to exact results indi- 
cate that UHF works well for ground state properties in 
random systems [TH 0|.A11 calculations presented were 
for sufficiently large systems (typically ~ 16 3 sites) that 
finite size effects are negligible. In Fig^ we show the 
single particle density of states and the distributions of 
both hole density and induced local moments for fixed 
U = 0.4 W,n h = 0.12, x = 4% and three values of 
V = 0.2, 0.45, 0.8 W. The Fermi energy is at zero. For 
small V there are no local magnetic moments. As V in- 
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FIG. 1: (Color online) Density of states (left) and distribu- 
tions of charge and magnetic moments (middle and right col- 
umn). The right column shows the densities projected onto 
oxygen neighboring vacancies (nn); the middle column other 
oxygens (not nn). The majority and minority bands(left col- 
umn) are continuous and dashed respectively. Energies are in 
units of bandwidth measured from the Fermi level, the mo- 
ments in units of fiB- 



creases, moments appear (V > 0.28 W) as the charge and 
spin distributions develop structures localised around the 
vacancies (right column). Later (V w 0.45 W), a well- 
defined impurity band for the minority band splits from 
the valence band. The first peak in the distributions 
(right column), at low densities, corresponds to sites with 
Vi = V/3 and the second broad peak to Vi = V. As V 
further increases, we also observe a splitting of the den- 
sity of states in the majority band and holes concentrated 
around the vacancies, providing saturated localized mo- 
ments. Note that the total moment on each cube centered 
on a vacancy is almost saturated (~ 3^_b) for V > 0.4 W. 
We now estimate the magnetic couplings at all distances 
between the different moments, by UHF calculation of the 
carrier Green functions Gf^ . We will subsequently make 
a self-consistent local random phase approximation (SC - 
LRPA) on the effective random Heisenberg model|l7|. 
which was successful for diluted magnetic semiconduc- 
tors. Each spin of the Heisenberg model will represent 
the total spin around one vacancy. Calculation of the 
exchanges is made by an extension of the method of 
Lichtenstein et a/.0]. They determined the exchange 
between any pair of (local) moments by calculating the 
change in energy when fields are applied to the sites of 
the two moments. Here, the fields are applied to two 
differing directions, the first to all the N$ = 8 neigh- 
bors i of one defect a; the second to the neighbors j of 
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FIG. 2: (Color online) Magnetic couplings (in Kelvin) as a 
function of separation of vacancies, in units of the cubic lattice 
spacing, for different values of the potential V. In each case 
U = 0AW, x = 0.04 and n h = 0.12 as in Fig.Q 



FIG. 3: (Color online) Curie temperature as a function of the 
potential V (in units of W). Other parameters as in Fig. 



the second b. Thus we calculate the exchange energies 
by summing the contribution of N$ x Ng pairs of terms 

Ja.b = Eieajeb where 



J, 



3 / ^Gl^i^Gl^du; 



3,1 ^ 



In UHF the local potential E< = U(nj -n\). 

The dilute Heisenberg model is denned with spins cen- 
tered on the position of the randomly placed vacancies 
and interacting by the J a ^ b . H = - J2 a ,b Ja,bS a • S b with 
S a the total spin around each vacancy. S a ~ 3[1b are 
large enough to be treated classically, allowing applica- 
tion of the local force theorem. In Fig|2]we show the cal- 
culated couplings averaged over different impurity con- 
figurations. We vary the potential V and fix x = 4% and 
rih/x = 3. This last choice (rather than 4) will be made 
clearer by the figures following. For small V the couplings 
oscillate with distance, but with antiferromagnetic near- 
est neighbor coupling. As V increases the couplings be- 
come more ferromagnetic but with further increase some 
become antiferromagnetic again. Thus from FigEI there 
is a range of values V where the couplings, while fluctuat- 
ing, are all ferromagnetic (J(r) > 0). This is associated 
in Figure^ (b) with the incipient development of a visible 
impurity band just at the band edge. This ferromagnetic 
bias in the "RKKY-like" oscillations corresponds to the 
resonant form of the impurity band 19] . 

We now treat the thermodynamics within SC- 
lrpa[I3. We recall that randomness is treated ex- 
actly, essential at the low concentrations relevant here. 
In Fig. we show the Curie temperatures (Tc) calcu- 
lated from the exchanges as a function of V. It is seen 
that there is a well-defined region where Tc is above 




FIG. 4: (Color online) Curie temperature for different carrier 
densities for 4% defects. V = 0AW, U = 0.4W are fixed. 



room temperature. Comparing to Figs. ^ and |3 we see 
that for V either too small or too large, Tc vanishes 
because the couplings are frustrated by RKKY-type os- 
cillations (small V) or superexchange (large V). As seen 
in Fig |3 near the optimal value (V=0.45W) the cou- 
plings are largest and do not change sign. Remark that 
at that value (see Fig ^b)) the impurity band is very 
similar to that of doped Ga(Mn)As|2fJ. Optimal values 
of V have high Tc for the same reason, the position of 
the impurity band, that Ga(Mn)As has a larger Tc than 
Ga(Mn)N and In(Mn)As which resemble more cases (c) 
and (a) respectively [lj . 

In Fig. 0| we fix the density of defects and V near the 
optimal value of Fig. [3] and vary the number of holes per 
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vacancy. Below and above critical values of jx there is 
no ferromagnetism. There is a window of concentrations 
where Tp becomes very large. The reason for instabil- 
ity outside the window is for low carrier densities, domi- 
nant antiferromagnetic superexchanges and for high den- 
sities, RKKY-like oscillations. Note that for the points 
where Tc = our calculation predicts that are still lo- 
cal moments, in agreement, for example, with ab-initio 
calculations [2j for defects in CaO, but no long-range fer- 
romagnetic order. Interestingly, the formal charges of 
vacancies in HTO2 or CaO are near the edges of stabil- 
ity of ferromagnetism for the parameters chosen. This 
may suggest an explanation for the extreme sensitivity 
of current results to sample history. Note that a more di- 
rect comparison with experiments would require a precise 
knowledge of the parameters and the concentrations of all 
defects including especially Oxygen vacancies. Oxygen 
vacancies, would decrease the effective hole density and 
could stabilize Tc for HTO2 . The figure also suggests that 
a promising avenue for stabilizing or increasing the Curie 
temperature experimentally may be to use non-magnetic 
substitution of Hf(Zr) by Li, Na, K, Rb or Cs .. rather 
than Hf(Zr) vacancies which are difficult to control. 

Our results are based on a Hartree-Fock approxima- 
tion that tends to overestimate the tendency to ferromag- 
netism in some homogeneous itinerant models. Nonethe- 
less KanamoriHII, Tasaki[I2 and others jHQ E3| have 
shown that in "flat band" models in which there is a 
peak in the density of states with a relatively flat back- 
ground, ferromagnetism is indeed possible. Because of 
the vacancies, our model resembles a fiat band model. 
We remark that degeneracy of orbitals is not required, 
as the moments are stabilized by the molecular fields of 
other moments. Note that screening effects which reduce 
the effective U by the t matrix |21| are important when 
U ^> W which is not the case here. In addition, we 
have verified tat the results vary little with t' jt and U in 
broad ranges, 0.2 < U/W < 0.6 and 0.1 < t'/t < 0.6. In 
contrast Tc is very sensitive to both V and n^/x which 
are the most relevant parameters. We have also seen 
that below typically 1% no ferromagnetism is possible. In 
ref. [l2j| the equilibrium density of cationic vacancies was 
estimated as 0.03%. Thus higher non equlilibrium con- 
centrations are very likely needed, as are quite possible 
in films. However a more promising avenue for "d ' fer- 
romagnetism would be the subtitutional approach. 

To conclude, we present a model which catches the es- 
sential ingredients for vacancy- or substitutional-induced 
ferromagnetism of the oxygen holes. Our calculations 
shows that Curie temperatures above room temperature 
are feasible for a few per cent of vacancies or substitu- 
tions. The sensitivity of "d°" samples may be explained 
by the proximity of the doping to the stability bound- 
ary. The underlying mechanism is close to that of diluted 
III-V semiconductors: enhanced ferromagnetic couplings 
between resonant impurity levels at parameters just at 



the point where the impurity states split off from the 
valence band. It would be interesting to measure the 
profile of the density of states (e.g. by photoemission) to 
see whether it does indeed correspond to the proposed 
picture (see Fig. Q] center row). As Curie temperatures 
are non-monotonic in both the strength of the potential 
and the doping, accurate characterization may be needed 
in order to find useful ferromagnetic materials. Ab-initio 
results for the effective Hamiltonian could help refine the 
choice of compound. Our calculation, Figure 0] suggests 
that for Hf0 2 (Zr0 2 ), site substitution of the Hf(Zr) by 
elements differing by 3 in formal charge, i.e. Group 1A of 
the periodic table would be a promising direction. This 
may be attractive since it can be realized in bulk as well 
as films. 
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